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Abstract. The projected (PME) and virial mass estimator (VME) are 
revisited and tested using an A^-body equilibrium system. It is found that 
the PME can overestimate the mass by « 40% if a cluster of galaxies is 
sampled only about its effective radius, -Re- The exact value of this error 
depends on the mass profile of the cluster. The VME can yield errors of 
~ 20% if we are also able to sample about the R^. Both estimators yield 
the correct total mass for an equilibrium system. Theoretically, if we 
know the total extent of the system, the VME provides an accurate mass 
at different radii, provided its gravitational potential term is correctly 
taken into consideration. The VME can provide acceptable masses for 
equilibrium systems with anisotropic velocity dispersions. 



1. Introduction 

Clusters of galaxies provide an important tool to study the large-scale structure 
of the universe (e.g. Peebles 1980), and several methods are used to determine 
their total mass and mass profile. Recently, some questions regarding the accu- 
racy and consistency of the PME and VME have been raised by several authors. 

Carlberg, Yee & Ellingson (1997, hereafter CYE) found that the VME over- 
estimates mass by about ~ 20% in clusters of galaxies, and attribute this dis- 
crepancy to the neglect of a surface pressure {'iPV) term in the continuous form 
of the virial theorem. This ?>PV term is starting to be used as a correction 
term when mass estimates from the VME are calculated (e.g. Girardi, et al. 
1998). Perhaps more critical is the situation in A^-body simulations where there 
is no ambiguity regarding positions and velocities of particles, and the obtained 
masses using either the PME or the VME do not agree with the true mass (e.g. 
Thomas &: Couchman 1992). 

Here we revisit the PME and the VME and test their consistency using 
an A^-body equilibrium system, and explain the discrepancies that arise in the 
relevant literature. A more comprehensive treatment of these matters, along 
with a determination of the mass profile of Coma cluster can be found elsewhere 
(Aceves & Perea 1998). We consider here, however, the performance of the VME 
for a system in equilibrium but with a non-isotropic velocity distribution. 
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Figure 1. Schematic representation of the use of the PME and the 

VME by apertures in a chister of galaxies; Rb is the physical boundary 
of the system, R is the radius of the aperture, and Rij is the inter- 
particle separation between two galaxies. 



2. Projected and Virial Mass Estimators 

For an isotropic and spherically symmetric gravitational system the PME and 
the VME are, respectively (e.g. Heisler et al. 1985; Perea, et al. 1990): 

where is the observed linc-of-sight velocity of a galaxy relative to the cluster 
mean, Ri is its projected radius from the center of the distribution, and Rij 
denotes the projected inter-galaxy separation. 

If we consider partial sampling of a system, a correction term for the PME 
has to be included (e.g. Haller & Melia 1996, Aceves & Perea 1998), namely: 

A ^ ^{v'^R) - M{r) = rV(r)a2(r) , (2) 

where M(r) is the true mass of the system, p{r) is the mass density and <Jr(r) 
the radial velocity dispersion. For the VME, CYE have suggested a surface 
pressure term 3PV = Anr^pa"^ as a correction to virial mass estimates to avoid 
overestimation (see also Girardi, et al. 1998). Wc will now test these estimators, 
and the suggested correction terms, for an A^-body equilibrium model. 

Throughout this work the units used are such that the total mass is M = 1, 
the effective radius is i?c = !> and G = 1. A Monte-Carlo realization of A = 10^ 
particles for a Dehnen 7 = model (Dehnen 1993) was used to test the mass 
estimators referred above. The effective radius in this model, i.e. the locus 
where the projected mass is one-half of the total mass, is R^ = 2.9036a , with 
a being a particular scale radius. We use here Dehnen model as an illustrative 
case, due mainly to the simplicity of its expression for the velocity dispersion, 
but we have considered other models and similar trends as those described here 
were found; although with different numerical values. 
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Figure 2. Results of the application of the PME and VME to an 
A^-body equilibrium system following a Dehnen (7 = 0) profile. 



In order to apply the PME and the VME we proceeded in a similar fashion 

as one would observationally do. We divided the complete system in different 
'apertures' of increasing radii, and applied both the PME and the VME just to 
the particles inside each of them; see Fig. 1. This approach mimics, particularly 
for the VME, the lack of information one may have on the total extent of a real 
cluster of galaxies, especially if large amounts of dark matter exist outside their 
optical boundary. 

In Fig. 2, our results arc shown. Figure 2a. shows the difference between 
Mp and M(r), and with respect to M{R); the theoretical value of A(r) is also 
displayed. In Fig. 2b. the differences of the VME applied by apertures, both in 
its projected and non-projected versions, with the true mass are shown. Almost 
perfect agreement with the theoretical expectations of Eq. (2) exists for the 
PME, except by the Poisson noise in the discrete A?^-body model. Prom the 
values of A(r) and the numerical experiments it follows that the PME may 
overestimates the true mass by ^ 20% if sampling is made only around « IRe, 
and if the mass profile is similar to a Dehnen model. For profiles similar to 
Hernquist (1990) or De Vaucouleurs (1948), i.e. 'cuspy', the error is smaller 
while for systems of a larger core, e.g. a King Modified (see Aceves & Perea 
1998), it can be up to ^ 40%. 

The discrepancy of the VME result applied by apertures with the true mass 
requires a somewhat different treat. CYE attribute this kind of error to a neglect 
of the 3PV term in the continuous virial theorem: 



/r rr 
r'^pal dr = -47rG / rpM{r) dr 



(3) 



where the surface pressure term is 3PF = Airr^pa^ . We can readily estimate 
the position of the maximum of the 3PV term for a Dehnen model, where 
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Figure 3. Mass estimates using the aperture VME, Ma, and that 
including the contribution to the potential energy from particles outside 
this aperture, Mq. The true mass of the Dehnen model, Md, is also 
indicated. Here, the non-projected version of the VME is used, but the 
same holds for the projected quantity. 

obtaining Tmax = (7 + a/65)o/8 ~ 0.65, yielding a maximum value of SPFmax ~ 
0.04 ; see Fig. 2b. This value is about the maximum overestimate in mass when 
the 3PV term is not taken into consideration in (3). 

Moreover, without the 3PV term the resulting mass from Eq. 3 would 
be MNSp(r) = 3ra^{r)/G, which overestimates the mass for r ^ R^, the near 
isothermal region of the model, and underestimates it for r ^ R^. For example, 
we have that at r = lOORc the mass is Mnsp ~ 0.6 instead of ~ 1 ; i.e. about a 
40% underestimate error is obtained for the total mass. 

The physical reason for the discrepancy in the aperture evaluation of the 
VME value with the true mass M{r) is the incorrect evaluation of the potential 
energy term, i.e. that involving J^^/^ij^ equation (1). The previous sum- 
mation has to be done for all j-th particles in the system, not only over those 
inside the aperture radius defined by the i-th particle (Aceves & Perea 1998). 
When this is done, we find very good agreement between the calculated and 
theoretical masses at different radii; see Fig. 3. This is basically a consequence 
of the long-range nature of the gravitational interaction. Hence, there is no 
need to introduce a 3PV term as a correction term when Eq. 1 for the VME 
is used, we only need to account for the mass outside the particular boundary 
we are considering to evaluate the mass. An upper bound to the error when no 
knowleadge of the 'real' extent of the cluster is also provided by (2). 

When the physical boundary of a cluster is unknown the correction factor is 
model dependent. Regarding real clusters, we are faced with a practical problem 
in astronomy. As we have seen, there is no problem for the VME to estimate 
the mass in an equilibrium system at different radii, provided that we actually 
know its total extent. But in astronomy, the physical boundary is a somewhat 
ambiguous matter, affected by the same sensitivity of our instruments and even 
by the criteria used to 'define' a cluster of galaxies. Nonetheless, one can use 
the difference in the mass estimated from the PME and from the VME to have 
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Figure 4. Cold collapse of an initially oc l/r distribution of mat- 
ter, and 2T/W = 0.1 . {Left) The initial configuration of the system. 
(Right) System's appearance after having reached a quasi-steady state, 
2T/W ~ 1, at i = 31.25 time units. The final system is ellipsoidal 
rather than spherical, and has a non-isotropic velocity distribution. 



an idea of the extent of the sampling done, as illustrated by Aceves & Perea 
(1998) for the Coma cluster of galaxies. This approach is more or less correct 
only when we can use other indicators, e.g. X-rays, to infer that the cluster is 
near equilibrium. Further, the asymptotic approach to a particular value of the 
integrated mass is also indicative that the system is near complete sampling. 

In iV-body simulations (e.g. Thomas & Couchman 1992) the discrepancy 
seems to arise due to the use of formulae for the VME and the PME suitable only 
for test particles; i.e. not valid for mutually interacting particles. Furthermore, 
the possible presence of anisotropies in the outer parts of the numerical clusters 
can affect the results especially those from the PME. 



3. An Anisotropic System 

In this section, as an illustrative example, we apply the virial theorem to an 
anisotropic system in equilibrium. The previous can give us an idea of the 
errors one might expect to commit under such conditions. For that purpose, we 
performed an AT-body simulation of a cold collapse (e.g. van Albada 1982). The 
initial system had a ratio of twice the kinetic to potential energy of 2T/W = 0.1, 
and the particles were distributed with a density law p oc l/r . We followed the 
evolution until the system reached a quasi-steady state, see Fig. 4. 

In Fig. 5 we show our results. On the left panel, the radial velocity dis- 
persion and the projected velocity dispersion are shown for the final snapshot of 
the simulation; the anisotropy function /3(r) = 1 — ag/a'^ is shown in the inset. 
Note how the system is close to being isotropic, /? ~ 0, at the inner regions while 
particles on radial orbits dominate the outer parts. 

We applied the mass estimators described in §2 to the final configuration 
which is in a quasi-steady state. The PME and the VME were applied assuming 
isotropy in the system, i.e. formulae in Eq (1) were used, to reflect in a sense the 
way one proceeds in estimating the mass in clusters where this assumption is 
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Figure 5. (Left) Radial velocity dispersion (filled squares), and pro- 
jected velocity dispersion (open squares) in the final system of the 
cold collapse in Fig. 4. The inset depicts the anisotropy parameter 
P{r) = 1 — (Tg/a,^ . (Right) Application of the different mass estimators 
considered here. Mvr corresponds to the non-projected form of the 
VME. 



usually made; albeit observations indicating that clusters are somewhat triaxial 
(e.g. Plionis, et al. 1991). For the virial theorem estimates, we considered the 
contribution to the potential energy from the particles outside a given radius, 
and considered its projected version, Mvme, and its non-projected one, Mvr- 
The true mass M(r) was calculated by counting the particles inside spheres of 
different radii in the A'^-body system. 

It is interesting to note that the non-projected virial theorem follows very 
close the true mass of the system, and the projected one overestimates it by 
10% throughout most of the system's extent; see right panel of Fig. 5. This 
further supports the thesis that the VME is not very sensitive to anisotropies 
in equilibrium systems. We are not considering here the possible role of a mass- 
spectrum or rotation in clusters. 

The virial theorem tends to be more sensitive to the equilibrium state of 
the system than to its anisotropy. Since Mvr oc 2T/W it follows that for "cold" 
systems {2T/W ^ 1) the mass will be underestimated, while for "hot" systems 
{2T/W ^ 1) an overestimate will result. Indeed, for our initial configuration 
the mass estimate provided by the virial theorem was Mvr(i = 0) ~ 0.1 . 

The PME is the most sensitive to anisotropies in the system, hence more 
propense to give unreliable mass estimates even for a quasi-equilibrium state. 
Therefore, the VME provides a better estimate of the mass at different radii than 
the PME in situations where galaxies have non-isotropic velocity dispersions, and 
no information on the distribution of orbits is known. 



4. Conclusions 
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1. The projected mass estimator (PME) and the virial mass estimator (VME) 
yield accurate results for the total mass of a gravitational system, provided 
that the whole system is in steady-state, spherically symmetric, and has 
an isotropic velocity distribution. 

2. When partial sampling is done, the error committed in the use of the PME 

and VME depends on the density profile, and their maxima occurs about 
the system's effective radius. For realistic profiles, the error can be up to 
^ 40% for the PME and ^ 20% for the VME. 

3. When the total extent of the system is known, the VME yields accurate 
mass estimates at different radii. This is achieved by accounting for the 
potential energy contributed from the particles outside a given radius. 

4. A surface pressure correction term cannot be used in the mass calculated 
from the VME, when partial sampling of a system is done, to obtain a 
correct mass. The correction term in this circumstance is dependent on 
the mass distribution of the cluster. 

5. The VME tends to yield a mass profile closer to the true value in systems 
with anisotropic velocity distributions, provided this is in quasi-steady 
equilibrium. 
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